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Abstract.  A  new  computational  model  based  on  Eu’s  generalized  hydrodynamics,  which  has  been  recently 
proposed  for  describing  the  motion  of  gases  in  non-equilibrium  state  and  is  shown  to  be  consistent  with 
the  second  law  of  thermodynamics,  is  presented.  The  general  understanding  of  Eu’s  generalized  hydrody¬ 
namics,  which  employs  the  cumulant  expansion  for  the  Boltzmann  collision  integral  instead  of  the  BGK 
approximation,  is  also  obtained  by  considering  three  fundamental  flows;  compressed  gas  in  shock  waves, 
expanding  gas,  and  velocity  shear  flow.  The  study  on  these  problems  reveals  that  Grad’s  equations  are 
similar  to  Eu’s  equations  in  the  slip  flow,  but  become  drastically  different  from  Eu’s  equations  in  shock 
structure  problem.  A  plausible  explanation  is  that  the  relaxation  time  approximation  may  be  insufficient  in 
modeling  the  extreme  nonlinearity  of  shock  structure  since  the  Boltzmann  collision  integral  plays  a  critical 
role  in  this  case.  Finally,  by  considering  the  microscale  channel  flow,  a  new  slip  boundary  condition  based 
on  Langmuir’s  theory  is  presented  that  predicts  a  trend  of  increasing  pressure  curve  nonlinearity  with  in¬ 
creasing  rarefaction  and  a  minimum  in  mass  flow  rate,  which  are  not  the  case  with  the  results  predicted  by 
the  first-order  Maxwell  slip  condition. 


INTRODUCTION 


Analysis  of  high  thermal  non-equilibrium  gas  flows  in  many  instances  requires  mathematical  models  be¬ 
yond  the  Navier-Stokes  equations.  Many  computational  models-either  fully  kinetic  (DSMC)  or  fluid  dynamic 
(moment  equations)-have  been  proposed,  but  it  turned  out  that  some  have  difficulty  in  solving  low-speed 
micro-scale  flow  and  others  suffer  non-trivial  problems:  violation  of  the  2nd  law  of  thermodynamics,  closure 
breakdown,  and  trouble  in  boundary  condition.  This  status  is  somewhat  surprising  since  the  fluid  dynamic 
approach  is  simply  to  solve  the  following  collision-free  hyperbolic  conservation  laws, 

^  J  U dV  +  j  F  •  ndS  =  0,  (1) 

where  S  represents  the  bounding  surface  of  the  control  volume  V,  and  U  represent  conserved  variables 
(p,  pu,  pE).  F  in  the  surface  integral  represents  the  flux  consisting  of  non-conserved  variables  (stress  and 
heat  flux)  and  conserved  variables.  Here  it  should  be  emphasized  that  these  laws  are  the  exact  consequence 
of  both  kinetic  theory  and  continuum  mechanics.  Only  after  some  approximations  are  made  to  non-conserved 
variables,  they  become  approximate  fluid  dynamic  equations,  for  example,  the  Navier-Stokes  equations.  It  was 
found,  however,  that  the  procedure  to  calculate  non-conserved  variables  is  not  simple,  and  indeed  it  is  generally 
believed  that  there  exist  no  high-order  equations  fully  consistent  with  the  macroscopic  and  phenomenological 
thermodynamics . 
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EU’S  GENERALIZED  HYDRODYNAMICS 


Towards  solving  this  dilemma,  an  alternative  and  yet  unknown  method  in  which  the  second  law  of  thermody¬ 
namics  is  satisfied  to  every  order  of  approximation,  the  Eu’s  generalized  hydrodynamics  (GH)  [1]  is  considered 
in  this  study  as  the  basis  of  a  new  computational  model  for  high  non-equilibrium  gas  flows.  The  essence  of  his 
theory  can  be  summarized  as  follows.  Starting  from  the  following  irreversible  thermodynamics 

TdV  =  dE  +  pdV  +  J2  ®d<  /p  >,  (2) 

which  is  equivalent  to  an  exponential  form  of  the  distribution  function 

/  =  /(o)e-Ext‘>>®h(“>-'‘]/*BT)  (3) 

a  thermodynamically  consistent  hydrodynamic  equations  can  be  developed.  $  in  the  equation  (2)  represents 
a  non-equilibrium  entropy  that  recovers  the  equilibrium  entropy  S  for  an  equilibrium  system.  represent 
high-order  moments  such  as  stress  and  heat  flux,  while  represent  tensors  whose  role  is  similar  in  the 
coefficients  of  conventional  expansion  methods.  A  leading  term  X  M  can  be  written  as  — II / ( 2p) ,  where  II 
represents  stress  tensor,  is  directly  related  to  through  the  relation  §(“)  =<  h^a'f  >.  p  is  the 
normalization  factor.  The  essence  of  the  Eu’s  generalized  hydrodynamics  can  be  found  in  an  exponential  form 
of  the  distribution  function  (3).  The  use  of  the  exponential  form  makes  it  possible  to  apply  the  cumulant 
expansion  and  ensures  that  the  approximated  distribution  functions  always  remain  non-negative.  Note  that 
/  can  never  be  negative  and  must  either  have  finite  bounds  in  velocity  space  or  tend  to  zero  as  a  particle 
velocity  tends  to  infinity.  This  point  can  be  highlighted  by  realizing  that  the  following  conventional  polynomial 
expansions  show  the  poor  convergence,  meaning  that  there  always  exists  a  region  in  which  the  function  can  be 
either  negative  or  unbounded. 

/  =  /<°)  +  Kn/W  +  Kn2/(2)  +  •  •  • .  (4) 

Thus,  it  is  the  exponential  form  of  the  distribution  function  that  distinguishes  the  Eu’s  GH  from  the  rest  of 
other  hydrodynamic  approaches. 

One  of  the  simple  form  of  GH  constitutive  relations  retaining  the  essence  of  theory  can  be  written  as  [2,3] 

ilq(cR)  =  nNS  +  [n  •  Vu]<2\  (5) 

Q  g(cR)  =  Qns  +  n  ■  qNs>  (6) 

where  IINS  =  -2t?[Vu](2),  Qns  =  -AVlnT,  U  =  Q  =  Vii  =  &  = 

TL:  n  +  Q  •  Q,  q(cR)  =  sin^^ .  Here  Q  represents  heat  flux  vector. 

A  constant  c  has  a  value  between  1.014  (Maxwellian)  and  1.223  (^=3),  where  u  is  the  exponent  of  the  inverse 
power  laws  of  gas  particles  and  s  =  l/2  +  2/(i/  —  1).  In  case  of  helium  gas,  c  becomes  1.046(17  =  14).  77  and  A  are 
the  Chapman-Enskog  viscosity  and  thermal  conductivity  .  They  can  be  expressed  as  q  =  Ts,  A  =  Ts+1.  For  a 
monatomic  gas,  7  =  5/3  and  Pr  =  2/3.  A  composite  number,  which  is  defined  as  Ng  =  r)u/(pL)  ~  Kn M  and 
appears  in  the  term  measures  the  magnitude  of  the  viscous  stress  relative  to  the  hydrostatic  pressure, 
so  that  it  indicates  the  level  of  departure  from  thermal  equilibrium.  When  the  new  constitutive  relations  (5) 
and  (6)  are  compared  with  Grad  [5]  and  Burnett  [7]  equations,  several  interesting  results  can  be  found.  In 
multi-dimensional  flows,  the  following  three  cases  can  be  considered  most  basic  flows. 


Shock  Wave  and  Gaseous  Expansion  When  the  gas  experiences  compression  or  expansion,  the  gradient 
of  normal  velocity  will  remain  either  negative  or  positive.  In  these  cases  (du/d x  only),  the  Eu’s,  Grad’s,  and 
Burnett’s  relations  for  the  normal  stress  reduce  to 


(1  +  nxajni 


n* 


4nT 


4-7H, 


nT 


(1  +  )n. 


(7) 


In  Fig.  1  (a),  Eu’s  and  Grad’s  results  for  the  normal  stress  are  plotted.  It  can  be  observed  that  Eu’s  relation 
is  well-defined  in  all  non-equilibrium  states.  Note  that  the  Grad’s  equation  has  a  non-removable  singularity 
at  n  =  4/7  in  positive  region  (shock  wave)  and  the  Burnett’s  equation  loses  one-to-one  correspondence  in 
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negative  region  (expanding  gas),  which  in  turn  implies  the  violation  of  positive  entropy  production  [7].  Even 
though  the  algebra  in  these  relations  is  quite  simple,  it  turns  out  that  the  results  are  very  general.  In  the 
previous  work  [4],  it  was  found  that  even  if  heat  flux  is  included,  the  general  characteristics  remain  unchanged. 
Furthermore,  an  independent  calculation  done  by  Karlin  et  al.  [8]  confirms  this  characteristics.  They  extended 
the  hydrodynamic  description  into  a  highly  non-equilibrium  domain  by  summing  all  the  relevant  subseries  of 
the  original  Chapmann-Enskog  expansion  from  the  Grad  equations  and  obtained  the  following  linear  relation 
for  Maxwellian  molecules, 


IIt.^.  — 


211, 


n. 


-  2  +  [(rUNS  -  2)2  +  8ft2  ]1/2 


(8) 


It  can  be  easily  shown  that  this  relation  is  very  similar  to  the  GH  relation,  even  though  some  details  are 
different. 

The  salient  feature  of  the  GH  equation  is  that  it  recovers  the  free-molecular  limit  in  gaseous  expansion  and 
shows  high  nonlinearity  in  gaseous  compression.  It  can  be  induced  from  this  nonlinearity  that  there  must 
exist  at  least  one  region  of  shock  wave  structure  in  which  the  actual  stress  is  greater  than  the  stress  calculated 
by  the  Navier-Stokes  equations,  irrespective  of  Mach  numbers.  It  is  exactly  this  nonlinearity  that  makes  the 
GH  solutions  for  shock  structure  in  strong  agreement  with  experimental  data  [6,2].  Notice  that  the  essential 
contribution  to  this  nonlinearity  comes  from  terms  q(cR)  and  [II  •  Vu](2),  which  represent  energy  dissipation 
by  particle  collision  and  the  strong  coupling  between  the  stress  and  velocity  gradient.  On  the  other  hand,  the 
work  by  Karlin  et  al.  implies  that  summation  of  all  the  relevant  terms  is  necessary  to  explain  the  extreme 
nonlinearity  of  shock  structure.  This  point  may  explain  a  consistent  finding  that  the  common  relaxation  time 
(BGK)  approximation  for  the  Boltzmann  collision  integral  yields  poor  results  in  high  Mach  shock  structure 
problem. 


Velocity  Shear  Flow  When  there  exists  the  gradient  of  shear  velocity  only  only),  the  Eu’s,  Grad’s 
and  Burnett’s  relations  reduce  to 


Hx3T  - 


2^LS/3 


2n2  /8 

-,  ftxa;  = - — ,  ftxx  =  -o.733n2 

q*(cR)  +  2Tllmj3'  1  +  2H  StfNS/3’ 

where  R2  =  3fta,a,(fta;a;  -  1)  and  with  the  same  stress  eclipse  in  Eu’s  and  Grad’s  cases, 

3 

ft*j/  =  Sign(ftj;yNS ) [  —  (1  +  fta:a:)fta:a;]  ^  • 


£J/NS> 


(9) 


It  can  be  shown  in  Fig.  1  (b)  that  the  behavior  of  stresses  becomes  very  different  from  the  Navier-Stokes 
description  in  highly  non-equilibrium  states.  The  normal  stress  is  generated  by  the  shear  velocity  gradient  and 
the  shear  stress  approaches  zero  as  the  velocity  gradient  increases.  This  asymptotic  behavior  means  that  the  gas 
slips  near  the  solid  wall.  As  a  result,  the  velocity-slip  phenomenon  can  be  explained  in  purely  hydrodynamical 
terms.  Notice  also  that  Eu’s  results  are  very  similar  to  Grad’s  results,  except  for  the  existence  of  a  nonlinear 
term  q(cR)  that  comes  directly  from  the  Boltzmann  collision  integral.  This  similarity,  which  is  not  the  case  in 
shock  structure  problem,  can  be  explained  by  the  existence  of  the  stress  eclipse.  It  must  be  noticed  that  the 
kinematic  relations  are  solely  responsible  for  this  stress  constraint,  not  the  collision  integral  terms.  Therefore 
it  is  not  surprising  to  find  out  that  the  BGK  approximation  yields  fairly  good  results  in  the  boundary  layer 
flow,  which  is  dominated  by  the  velocity  shear  flow. 


Multidimensional  Flow  Even  though  the  new  equations  can  explain  nonlinear  transport  of  gas  in  non¬ 
equilibrium  state,  a  proper  computational  method  must  be  developed  in  order  that  they  may  serve  as  a  tool 
for  the  simulation  of  gas  flows  in  multidimensional  geometry.  In  the  previous  work  [3],  an  idea  to  overcome  the 
difficulty  of  multidimensional  extension  caused  primarily  by  the  existence  of  highly  nonlinear  terms  has  been 
proposed.  On  the  basis  of  such  idea,  a  numerical  solution  of  the  equations  (1),  (5),  and  (6)  was  obtained  for 
the  two-dimensional  flow  over  a  flat  plate  with  M  =  12.9  and  Kn  =  0.0067.  The  important  finding  was  that 
the  GH  calculation  produces  results  very  close  to  physical  intuition,  removing  the  singularity  near  the  leading 
edge  that  is  ill-defined  in  continuum  gas  dynamics.  The  ultimate  reason  for  this  can  be  traced  to  the  fact  that 
the  shear  stress  measured  by  the  GH  equations,  shown  in  Fig.  1  (b),  are  nonlinearly  related  to  the  gradient  of 
velocity,  smaller  than  that  measured  by  the  NS  equations  in  high  non-equilibrium  flow  regions. 
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FIGURE  1.  (a)  Eu’s  GH  relation  relative  to  the  Navier-Stokes  relation  in  the  region  of  compressed  (positive)  and 
expanding  (negative)  gas  (no  heat  flux  for  simplicity).  The  horizontal  axis  represents  the  velocity  gradient  (top),  (b) 
The  normal  and  shear  stresses  described  by  Eu’s  GH  and  Grad’s  equations  (bottom). 


MEASURE  OF  THERMAL  NON-EQUILIBRIUM 

According  to  the  conventional  non-equilibrium  theory  based  on  the  mean  free  path,  gas  flows  encountered  in 
low-density  environments  and  microscale  geometries  at  standard  atmospheric  conditions  can  be  treated  with 
no  difference  and  requires  the  mathematical  models  beyond  the  continuum  description.  Therefore,  it  is  natural 
to  argue  that  DSMC  or  high-order  moments  equations  should  be  used  to  simulate  microscale  gas  flows  with 
high  Knudsen  numbers.  There  exists,  however,  an  ambiguity  in  using  the  mean  free  path  as  the  parameter 
to  classify  non-equilibrium  gas  flow  regimes.  For  example,  there  are  many  calculations  in  which  the  problem 
with  a  smaller  Knudsen  number  shows  high  non-equilibrium  effects  than  the  problem  with  a  higher  Knudsen 
number.  In  fact,  the  classification  based  on  the  Knudsen  number  is  largely  empirical  and  depends  highly  on 
the  problem  considered.  This  is  not  surprising  at  all  if  we  recall  the  physical  reasoning  that  the  gas  flows 
associated  with  microscale  geometry  at  normal  pressure,  which  is  typical  in  most  of  MEMS  devices,  is  perfectly 
near  equilibrium  unless  the  gas-surface  interaction  has  long-range  effects  on  the  overall  flow  structure.  It  is 
well  known,  however,  that  the  scale  length  associated  with  the  gas-surface  interaction  (due  to  physical  force) 
is  order  of  molecular  size  (nanoscale). 

This  reasoning  may  give  a  hint  of  a  possibility  of  developing  a  better  measure  of  thermal  non-equilibrium 
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in  gas  flows.  Indeed,  if  the  GH  equations  (5)  and  (6)  are  carefully  examined,  it  becomes  immediately  evident 
that  a  composite  number  Ns,  which  is  defined  as 


Ns  - 


TjrUr 

prL 


' yM 2 

Re 


=  KnM 


indicates  the  level  of  departure  from  thermal  equilibrium.  When  Ns  is  small,  the  GH  constitute  relations 
recover  the  Navier-Stokes  relations.  This  can  also  be  confirmed  from  the  graphical  representation  of  the  GH 
relations  (Fig.  1). 

Even  though  this  finding  appears  trivial,  it  has  huge  implications.  The  most  obvious  one  is  that  the  classi¬ 
fication  of  bulk  flow  regimes  should  be  based  on  the  Knudsen  number  times  the  Mach  number,  not  Knudsen 
number  alone.  In  order  to  explain  this,  let  us  consider  a  typical  rarefied  flow,  for  example,  a  hypersonic  flow 
over  a  vehicle  flying  at  an  altitude  of  93  km.  For  realistic  values  of  flight  conditions  (length  =  lm,  velocity= 
8 km/s),  dimensionless  parameters  can  be  approximated  as 


Kn  =  0.037,  Ns  =  1.0. 


On  the  other  hand,  for  helium  gas  flow  in  a  typical  microchannel  (length  =  1.2 pm,  exit  velocity  =  50cm, /s), 
they  are 


Kn  =  0.145,  Ns  =  7.67  x  10“5. 

According  to  the  new  picture,  there  is  no  need  of  using  the  high-order  equations  in  microscale  gas  flows  since 
Ns  is  extremely  small.  Notice  that  Ns  will  remain  relatively  small  for  low-speed  flow  even  when  the  channel 
height  is  measured  on  a  nanoscale.  The  Navier-Stokes  relations  will  be  just  enough  to  calculate  these  flows  if 
the  slip  effect  on  the  solid  surface  due  to  high  Knudsen  number  is  properly  taken  into  account.  However,  this 
line  of  thinking  has  not  been  validated,  primarily  due  to  the  lack  of  proper  slip  boundary  conditions. 


MICROSCALE  GAS  FLOWS 

The  lack  of  proper  slip  boundary  condition  can  be  manifested  when  a  benchmark  problem  is  considered, 
namely,  internal  gas  flows  in  microchannels.  For  the  sake  of  simplicity,  let  us  consider  a  long  microchannel 
with  the  ratio  of  the  channel  height  ( H )  to  its  length  (L)  e  in  isothermal  condition  (273 K).  Using  the  method 
described  by  Arkilic  et  al.  [11],  the  Navier-Stokes  equations,  at  the  zeroth  order,  may  be  written  as 

_9(pu)  d(pv)  _ 
dx  dy 

e  dp  d2u 
Ns  dx  dy2 

The  streamwise  coordinate  x  and  the  wall-normal  coordinate  y  are  nondimensionalized  by  L  and  H,  respectively. 
The  reference  state  is  chosen  as  the  outlet  conditions.  The  reference  velocity  is  the  area-average  streamwise 
velocity  at  the  channel  exit.  The  solutions  to  these  equations  can  be  obtained  by  either  numerical  methods  or 
analytical  methods  if  boundary  conditions  on  the  solid  wall  are  specified. 

Maxwell  Slip  Boundary  Conditions  The  Maxwell  slip  boundary  condition  has  been  popular  in  the 
past  owing  to  its  simplicity.  It  is  largely  based  on  the  notion  of  accommodation  coefficients  which  measure  the 
fraction  of  molecules  undergoing  diffusive  reflection.  It  turned  out  in  the  present  problem  that  the  first-order 
Maxwell  slip  condition, 
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yields  simple  analytic  solutions  of  velocity  profile  and  mass  flow  rate, 
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Here  <5  =  24 Ng/e,  px= o  =  Pin,  Px= 1  =  1,  and  a  represents  the  streamwise  momentum  accommodation  and 
is  often  taken  as  a  =  1.  It  was  found  from  these  solutions  that  the  prediction  by  the  first-order  Maxwell 
condition  is  in  strong  agreement  with  experimental  observations,  in  particular,  the  streamwise  velocity  profile 
and  the  nonlinear  pressure  distribution  along  the  channel.  However,  it  was  also  found  that  this  condition  cannot 
predict  the  existence  of  a  minimum  in  mass  flow  rate  in  the  transition  region,  which  highlights  the  essential 
nonlinearity  of  slip  phenomenon.  In  addition,  a  careful  examination  revealed  qualitatively  different  results 
with  experimental  findings  and  physical  intuitions.  For  example,  it  predicted  a  pressure  distribution  with 
smaller  nonlinearity  in  more  rarefied  flows,  suggesting  that  rarefaction  negates  compressibility,  in  contradiction 
to  experimental  results  and  physical  intuition.  It  also  predicted  a  drift  of  mass  toward  the  wall  as  the  flow 
progresses  down  the  channel,  in  contrast  with  the  physical  picture  which  can  be  induced  from  approximating 
the  channel  flow  by  Blasius  profiles  along  two  flat  plates. 

In  order  to  overcome  these  limitations,  the  high-order  Maxwell  conditions  have  been  proposed  by  Beskok  et 
al.  [13].  It  was  shown  that  the  high-order  conditions  predict  the  existence  of  a  minimum  in  mass  flow  rate.  It 
is  not  clear,  however,  that  they  can  predict  the  free-molecular  limit  of  mass  flow  rate  and  solve  the  problems 
described  in  previous  paragraph. 


Langmuir  Slip  Boundary  Condition  Even  though  the  description  of  slip  phenomenon  by  Maxwell  is 
simple,  there  exist  non-trivial  issues.  For  example,  it  lacks  the  predictability  since  this  condition  requires  some 
adjustable  parameters  in  the  form  of  accommodation  coefficients,  which  is  highly  dependent  of  type  of  gases  and 
wall  materials  and  geometry.  It  can  also  be  noticed  that  the  slip  velocity  by  this  condition  is  unbounded,  which 
may  cause  unphysical  high  slip  velocity  and  affect  the  type  of  vorticity  near  the  wall.  Finally,  it  is  susceptible 
to  the  numerical  errors  by  under-resolved  calculation  of  the  derivatives  of  tangential  velocity  near  the  wall. 
Here  a  simple  slip  boundary  condition  based  on  the  theory  of  adsorption  phenomena  developed  by  Langmuir 
[9,10]  will  be  considered.  This  condition  has  been  implemented  in  the  simulation  of  rarefied  hypersonic  flow 
over  a  flat  plate  [3]  and  microscale  gas  flows  [4]. 

This  condition  is  based  on  the  theory  of  adsorption  phenomena  pioneered  by  Langmuir.  In  his  historic  work 
of  molecular  films  of  gases  on  solid  surfaces,  he  developed  a  conception  of  adsorption  phenomena.  According  to 
his  description,  gas  molecules  do  not  in  general  rebound  elastically,  but  condense  on  the  surface,  being  held  by 
the  field  of  force  of  the  surface  atoms.  These  molecules  may  subsequently  evaporate  from  the  surface,  resulting 
in  some  time  lag.  Adsorption  and  slip  is  the  direct  result  of  this  time  lag.  From  this  physical  reasoning, 
he  derived  the  following  equation  of  the  fraction  of  surface  covered  at  equilibrium,  which  is  of  far-reaching 
significance  in  surface  chemistry, 
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Here  an  expression  of  /3, 
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can  be  introduced.  A  is  the  mean  area  of  a  site,  De  is  the  potential  parameter  and  they  can  be  inferred  from 
either  theoretical  prediction  or  experimental  data.  I  is  the  mean  free  path,  and  kB  is  the  Boltzmann  constant. 
£  is  a  mean  collision  distance  between  the  wall  surface  and  the  gas  molecules  at  all  angles.  With  a  calculated, 
the  boundary  values  of  dimensionless  velocity  can  be  determined  by 
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in  the  present  case, 


where  0  =  1.28/Kn  for  He-Al. 

It  can  be  easily  shown  that  this  new  condition  yields  the  following  results  [4], 
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satisfies  a  nonlinear  equation 
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Surprisingly,  it  can  be  shown  that  these  solutions  explain  all  the  contradictory  problems.  The  new  theory 
predicts  that  the  nonlinearity  in  pressure  distribution  increases  as  the  Knudsen  number  increases.  In  addition, 
the  new  theory  shows  that  there  always  exists  a  drift  of  mass  from  the  wall  and  its  magnitude  is  much  smaller 
than  one  predicted  by  the  Maxwell  slip  condition.  Lastly,  it  can  be  easily  shown  in  the  relation  (16)  that  there 
always  exists  a  minimum  in  mass  flow  rate  in  transition  regime.  This  critical  property  can  be  seen  in  Fig. 
2  where  the  flow  rate  exhibits  a  correct  continuum  limit  and  an  algebraic  singularity  at  very  high  Knudsen 
number, 
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In  fact,  the  mass  flow  rate  can  be  expressed  as  three  parts;  continuum  and  slip  (slip  +  self-diffusion). 

rh  —  nicom  inuum  T  fllglip  m8elf— diffusion > 


(17) 


(18) 


where 


^slip  —  j 


2  r  pL  - 1 


dp 

+  di 


x=1j 


^■self— diffusion 


1 

i  +  p' 


Note  that  m  self-diffusion  is  independent  of  the  pressure  difference  and  remains  finite  in  non-continuum  regime 
even  for  constant  pressure.  This  conclusion  agrees  well  with  the  explanation  advanced  by  Liu  [12]. 


SUMMARY 

A  new  computational  model  fully  consistent  with  the  second  law  of  thermodynamics  is  presented.  Through 
a  critical  review  of  the  previous  models  on  basic  flows,  the  characteristics  of  the  new  model  are  identified. 
The  important  finding  is  that  the  level  of  thermal  non-equilibrium  can  be  measured  by  the  composite  number 
Ns  and  consequently  the  Navier-Stokes  equations  with  proper  slip  boundary  conditions  can  describe  most  of 
microscale  gas  flows.  This  point  is  demonstrated  by  developing  the  new  slip  boundary  condition  based  on  the 
Langmuir  theory  of  the  adsorption  of  gases  on  metals  as  monatomic  films.  It  is  shown  that  all  the  critical 
features  including  the  existence  of  the  minimum  in  mass  flow  rate  can  be  explained  by  this  simple  condition. 

It  is  also  worthwhile  mentioning  that  the  approach  taken  in  the  present  study  holds  for  diatomic  gases  and 
dense  gases  and  therefore  the  essence  of  the  present  work  can  go  over  to  these  gases.  The  results  to  these 
problems  will  be  reported  in  due  course. 


ACKNOWLEDGMENTS 

This  work  was  performed  while  the  author  held  a  visiting  scholarship  in  Institute  for  Mathematics  and 
its  Applications  of  University  of  Minnesota  (1999-2000  Program:  Reactive  Flow  and  Transport  Phenomena). 
The  work  was  also  supported  in  part  by  Korea  Science  and  Engineering  Foundation  under  Research  Grant 
1999-2-305-001-3.  The  author  is  thankful  to  Professor  B.  C.  Eu  for  the  discussion  of  his  GH  equations. 


REFERENCES 

1.  Eu,  B.  C.,  Kinetic  Theory  and  Irreversible  Thermodynamics ,  Wiley,  New  York,  1992. 

2.  Myong,  R.  S.,  Phys.  Fluids  11-9,  2788  (1999). 

3.  Myong,  R.  S.,  AIAA  Paper  99-3578  (1999). 

4.  Myong,  R.  S.,  AIAA  Paper  2000-2500  (2000). 

5.  Grad,  H.,  Comm.  Pure  Appl.  Math.  2,  331  (1949). 

6.  Al-Ghoul,  M.,  and  Eu,  B.  C.,  Phys.  Rev.  E  56-3,  2981  (1997). 

7.  Comeaux,  K.  A.,  Chapman,  D.  R.,  and  MacCormack,  R.  W.,  AIAA  Paper  95-0415  (1995). 

8.  Karlin,  I.  V.,  Dukek,  G.,  and  Nonnenmacher,  T.  F.,  Phys.  Rev.  E  55-2,  1573  (1997). 

9.  Langmuir,  I.,  J.  Am.  Chem.  Soc.  40,  1361  (1918). 

10.  Langmuir,  I.,  Chem.  Revs.  13,  147  (1933). 


311 


FIGURE  2.  (a)  Variation  of  mass  flow  rate  as  a  function  of  exit  Knudsen  number  in  the  microchannel  (top),  (b) 
Variation  of  components  (continuum  and  slip)  of  mass  flow  rate  in  Langmuir’s  boundary  condition  (bottom). 


11.  Arkilic,  E.  B.,  Schmidt,  M.  A.,  and  Breuer,  K.  S.,  J.  Microelectromechanical  Systems  6-2,  167  (1997). 

12.  Liu,  C.  Y.,  Phys.  Fluids  11-3,  481  (1968). 

13.  Beskok,  A.,  and  Karniadakis,  G.  M.,  Microscale  Therm.  Eng.  3-1,  43  (1999). 


312 


